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I. INTRODUCTION 

Orbitally excited states of the D meson, denoted here 
as Dj, where J is the spin of the meson, provide a unique 
opportunity to test the Heavy Quark Effective Theory 
(HQET) [1, 2]. The simplest Dj meson consists of a 
charm quark and a light anti-quark in an orbital angular 
momentum L = 1 (P-wave) state. Four such states are 
expected with spin-parity J p = + (j = 1/2), 1+ (j = 
1/2), 1+ (j = 3/2) and 2+ (j = 3/2), which arc labeled 
here as Dq, D[, Di and D2, respectively, where j is a 
quantum number corresponding to the sum of the light 
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quark spin and the orbital angular momentum L. 

The conservation of parity and angular momentum in 
strong interactions imposes constraints on the strong de- 
cays of Dj states to Dtt and D*ir. The j = 1/2 states 
are predicted to decay exclusively through an S-wave: 
£>o -> Dtt and D[ -> D*n. The j = 3/2 states are 
expected to decay through a D-wavc: £>i — > D*7T and 
— > _D7r and Z?*7r. These transitions are summarized 
in Fig. 1. Because of the finite c-quark mass, the two 
J p = 1 + states may be mixtures of the j = 1/2 and 
j = 3/2 states. Thus the broad D[ state may decay 
via a D-wave and the narrow D\ state may decay via an 
S-wave. The j = 1/2 states with L = 1, which decay 
through an S-wave, are expected to be wide (hundreds of 
MeV/c 2 ), while the j — 3/2 states that decay through a 
D-wave are expected to be narrow (tens of MeV/c 2 ) [2- 
4]. Properties of the L = 1 Dj mesons [5] are given in 
Table I. 

The narrow Dj mesons have been previously observed 
and studied by a number of experiments [6-16]. Dj 
mesons have also been studied in semileptonic B de- 
cays [17-24]. Precise knowledge of the properties of the 
Dj mesons is important to reduce uncertainties in the 
measurements of semileptonic decays, and thus the de- 
termination of the Cabibbo-Kobayashi-Maskawa [25] ma- 
trix elements \V c b\ and |V^f,|. The Belle Collaboration 
has reported the first observation of the broad Dq° and 
Df mesons in B decay [12]. The FOCUS Collaboration 
has found evidence for broad structures in D + tt~ final 
states [13] with mass and width in agreement with the 
Dq° found by Belle Collaboration. However, the Parti- 
cle Data Group [5] considers that the J and P quantum 
numbers of the Dq° and D'® states still need confirma- 
tion. 

In this analysis, we fully reconstruct the decays B~ — > 
Z) + 7r~7r~ [26] and measure their branching fraction. We 
also perform an analysis of the Dalitz plot (DP) to mea- 
sure the exclusive branching fractions of B~ — ■* D°jTt~ 
and study the properties of the Dj mesons. The decay 
B~ — > D + tt~tt~ is expected to be dominated by the in- 



L=0- 



L= 1 



termediate states D^tt 



and D 



and has a possible 



contribution from B~ — > D + n~7t~ nonresonant (NR) de- 
cay. The D® and D'± states can not decay strongly into 
Dtt because of parity and angular momentum conserva- 
tion. However, the £>*(2007)° (labeled as D* here) mass 
is close to the Dn production threshold and it may con- 
tribute as a virtual intermediate state. The B* (labeled 
as B* here) produced in a virtual process B~ — * B*tt~ 
may also contribute via the decay B* — > D + ir~ . Possible 
contributions from these virtual states are also studied in 
this analysis. 
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FIG. 1: Mass spectrum for cu states. The vertical bars show 
the widths. Masses and widths are from Ref. [5]. The dotted 
and dashed lines between the levels show the dominant pion 
transitions. Although it is not indicated in the figure, the two 
1 + states may be mixtures of j — 1/2 and j — 3/2, and D[ 
may decay via a D-wave and Di may decay via an S-wave. 



TABLE I: Properties of L = 1 D°j mesons [5]. 





J p 


Mass Width 


Decays 


Partial 






(MeV/c 2 ) (MeV) 


seen [5] 


waves 


u a 


0+ 


2352 ± 50 261 ± 50 


Dtt 


S 


Df 


1+ 


2427 ± 36 384+^° 


D*n 


S, D 


Dl 


1+ 


2422.3 ± 1.3 20.4 ±1.7 D* 


7T, D°7T + 7r" 


S, D 


Df 


2+ 


2461.1 ± 1.6 43 ±4 


D*tt, Dtt 


D 



pie consists of 347.2 ffcT 1 corresponding to (382.9±4.2) x 
10 6 BB pairs (N B g) taken on the peak of the T(4S) res- 
onance. Monte Garlo (MC) simulation is used to study 
the detector response, its acceptance, background, and 
to validate the analysis. We use GEANT4 [27] to sim- 
ulate resonant e + e~ — > T(4S) — > BB events (generated 
by EvtGen [28]) and e + e _ — > qq (where q = u, d, s or c) 
continuum events (generated by JETSET [29]). 

A detailed description of the BABAR detector is given 
in Ref. [30]. Charged particle trajectories are measured 
by a five-layer, double-sided silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) immersed in a 1.5 T 
magnetic field. Charged particle identification (PID) is 
achieved by combining information from a ring-imaging 
Cherenkov device with ionization energy loss (dE/dx) 
measurements in the DCH and SVT. 



III. EVENT SELECTION 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
storage rings at SLAC between 1999 and 2006. The sam- 



Five charged particles are selected to reconstruct de- 
cays of B~ -> D+n-n- with D+ -> K-tt+tt+. The 
charged particle candidates are required to have trans- 
verse momenta above lOOMeV/c and at least twelve hits 
in the DCH. A K~ candidate must be identified as a 
kaon using a likelihood-based particle identification algo- 
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FIG. 2: K + -k~-k~ invariant mass distribution for D + candi- 
dates for the selected B~ —> D + it~-k~ decays without the 
cut on the mass of D + . Data (points with statistical errors) 
are compared to the results of the fit (solid curve), with the 
background distribution marked as a dashed line. The shaded 
area marks the D + signal region. 

rithm (with an average efficiency of ^85% and an average 
misidentification probability of ^3%). Any combination 
of K~tt + tt + candidates with a common vertex and an 
invariant mass between 1.8625 and 1.8745 GeV/c 2 is ac- 
cepted as a D + candidate. We fit the invariant mass 
distribution of the K ~tt + tt + candidates with a function 
that includes a Gaussian component for the signal and a 
linear term for the background. The signal parameters 
(mean and width of Gaussian) and slope of the back- 
ground function are free parameters of the fit. The data 
and the result of the fit are shown in Fig. 2. The invariant 
mass resolution for this D + decay is about 5.2 MeV/c 2 . 
The B~ candidates are reconstructed by combining a D + 
candidate and two charged tracks. The trajectories of 
the three daughters of the B~ meson candidate are con- 
strained to originate from a common decay vertex. The 
D + and B~ vertex fits are required to have converged. 

At the T(4S) resonance, B mesons can be character- 
ized by two nearly independent kinematic variables, the 
beam-energy substituted mass Toes and the energy dif- 
ference AE: 

™es = ^{s/2+p a -p B f/El-pl, (1) 
AE = E* B - (2) 

where E and p are energy and momentum, the subscripts 
and B refer to the e + e~-beam system and the B can- 
didate, respectively; s is the square of the center-of-mass 
energy and the asterisk labels the center-of-mass frame. 
For B~ — > D + Tt~Tr~ signal events, the toes distribution 
is well described by a Gaussian resolution function with 
a width of 2.6 MeV/c 2 centered at the B~ meson mass, 
while the AE distribution can be represented by a sum of 
two Gaussian functions with a common mean near zero 
and different widths with a combined RMS of 20 MeV. 
Continuum events are the dominant background. Sup- 



pression of background from continuum events is pro- 
vided by two topological requirements. In particular, we 
employ restrictions on the magnitude of the cosine of the 
thrust angle, cos Oth, defined as the angle between the 
thrust axis of the selected B candidate and the thrust 
axis of the remaining tracks and neutral clusters in the 
event. The distribution of | cos 0th | is strongly peaked to- 
wards unity for continuum background but is uniform for 
signal events. We also select on the ratio of the second to 
the zeroth Fox- Wolfram moment [31], R2, to further re- 
duce the continuum background. The value of R2 ranges 
from to 1. Small values of R2 indicate a more spher- 
ical event shape (typical for a BB event) while values 
close to 1 indicate a 2-jet event topology (typical for a 
qq event). We accept the events with | cos 0^1 < 0.85 
and i?2 < 0.30. The |cos0 t ^| (R2) cut eliminates about 
68% (71%) of the continuum background while retaining 
about 90% (83%) of signal events. 

To suppress backgrounds, restrictions are placed on 
to ES : 5.2754 < to es < 5.2820 GeV/c 2 , and AE: -130 < 
AE < 130 MeV. The selected samples of B candidates 
are used as input to an unbinned extended maximum 
likelihood fit to the AE distribution. The result of the 
fit is used to determine the fractions of signal and back- 
ground events in the selected data sample. For events 
with multiple candidates (~ 3.5% of the selected events) 
satisfying the selection criteria, we choose the one with 
best x 2 from the B vertex fit. Based on MC simulation, 
we determine that the correct candidate is selected at 
least 65% of the time. We fit the Toes distribution of 
the selected B~ — ► D + ir~Tr~ candidates with a sum of a 
Gaussian function for the signal and a background func- 
tion for the background having the probability density, 
P(x) oc xyl — x 2 exp(— £(1 — a; 2 )), where x — toes/to 
with too fixed at 5.29 GeV/c 2 and £ is a shape param- 
eter [32]. The signal parameters (mean, width of Gaus- 
sian) and the shape parameter of the background func- 
tion are free parameters of the fit. The data and the 
result of the fit are shown in Fig. 3a. We fit the AE dis- 
tribution of the selected B~ — ► D + ir~Tr~ candidates with 
a sum of two Gaussian functions with a common mean 
for the signal and a linear function for the background. 
The signal parameters (mean, width of wide Gaussian, 
width and fraction of narrow Gaussian) and the slope of 
the background function are free parameters of the fit. 
The data and the result of the fit are shown in Fig. 3b. 
The resulting signal yield is 3496 ± 74 events, where the 
error is statistical only. A clear signal is evident in both 
Toes and AE distributions. 

To distinguish signal and background in the Dalitz plot 
studies, we divide the candidates into three subsamples: 
the signal region, -21 < AE < 15 MeV, the left side- 
band, -109 < AE < -73 MeV, and the right sideband, 
67 < AE < 103 MeV. The events in the signal region are 
used in the Dalitz plot analysis, while the events in the 
sideband regions are used to study the background. 

In order to check the shape of the background AE dis- 
tribution, we have generated a background MC sample 
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FIG. 3: (a) jties and (b) AE distributions for D + ir~iv~ can- 
didates. Data (points with statistical errors) are compared 
to the results of the fits (solid curves), with the background 
contributions marked as dashed lines. The histograms are the 
corresponding distributions of the background MC sample as 
described in the text. The shaded area in (a) shows the sig- 
nal region, while the three shaded areas in (b) mark the signal 
region in the center and the two sidebands. 



of resonant and continuum events with B~ 



7T 7T 



signal events removed. The background MC sample has 
been scaled to the same luminosity as the data. The AE 
distribution of the selected events from the background 
MC sample is shown as the histogram in Fig. 3b. A 
small amount of peaking background is found from mis- 
reconstructed decays of B — > D + p~ with p~ — ► 7r~7r°, 
where a ir° is missed and a random track in the event is 
misidentified as a signal 7r _ . The background histogram 
in Fig. 3b is fitted with a sum of two Gaussian func- 
tions with a common mean for the peaking background, 
with parameters fixed to those obtained from the fit to 
data, and a linear function to describe the combinatorial 
background. The amount of peaking background is esti- 
mated at 82 ± 41 events. After peaking background sub- 
traction, the number of signal events above background 
is iV s ig na i = 3414 ± 85. The background fraction in the 
signal region is (30.4 ± 1.1)%. 



IV. DALITZ PLOT ANALYSIS 

We refit the D + and B~ candidate momenta by con- 
straining the trajectories of the three daughters of the B~ 
meson candidate to originate from a common decay ver- 
tex while constraining the invariant masses of K~tt + tt + 
and D + ir~ir~ to the D + and B~ masses [5], respectively. 
The mass-constraints ensure that all events fall within 
the Dalitz plot boundary. 

In the decay of a into a final state composed 
of three pseudo-scalar particles (D + n~n~), two degrees 
of freedom are required to describe the decay kinemat- 
ics. In this analysis we choose the two Dtt invariant 
mass-squared combinations x — m 2 (D + ir^) and y — 
to 2 (_D + 7t^~) as the independent variables, where the two 
like-sign pions tt^ and ir^ are randomly assigned to x and 
y. This has no effect on our analysis since the likelihood 
function (described below) is explicitly symmetrized with 
respect to interchange of the two identical particles. 

The differential decay rate is generally given in terms 
of the Lorentz-invariant matrix element M by 



dxdy 



\M\< 



2567r 3 m| ' 



(3) 



where tub is the B meson mass. The Dalitz plot gives 
a graphical representation of the variation of the square 
of the matrix element, |.M| 2 , over the kinematically ac- 
cessible phase space (x,y) of the process. Non-uniformity 
in the Dalitz plot can indicate presence of intermediate 
resonances, and their masses and spin quantum numbers 
can be determined. 



A. Probability Density Function 

We describe the distribution of candidate events in the 
Dalitz plot in terms of a probability density function 
(PDF). The PDF is the sum of signal and background 
components and has the form: 



PDF(x,y) 



/bg 

(1 - /bg 



B(x,y) 



fopB(x,y)dxdy 

[S(x,y)®K]e(x,y) 

/dp \- s ( x > v) ® n \ e ( x > y) dxd v ' 

(4) 



where the integral is performed over the whole Dalitz 
plot, the S(x,y) CEP 72. is the signal term convolved with 
the signal resolution function, B(x,y) is the background 
term, /b g is the fraction of background events, and e is 
the reconstruction efficiency. 

An unbinned maximum likelihood fit to the Dalitz plot 
is performed in order to maximize the value of 



C j f PDF(a i)W ) 



(5) 
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with respect to the parameters used to describe S, where 
Xi and yi are the values of x and y for event i respectively, 
and -/V event is the number of events in the Dalitz plot. In 
practice, the negative-log-likelihood (NLL) value 

NLL = -ln£ (6) 

is minimized in the fit. 

B. Goodness-of-fit 

It is difficult to find a proper binning at the kinematic 
boundaries in the x-y plane of the Dalitz plot. For this 
reason, we choose to estimate the goodness-of-fit x 2 i n 
the cos 9 (range from -1 to 1) and m^ lin (£)7r) (range from 
4.04 to 15.23 GcV 2 /c 4 ) plane, which is a rectangular rep- 
resentation of the Dalitz plot. The parameter 9 is the 
helicity angle of the Dtt system and Tn^^Dn) is the 
lesser of x and y. The helicity angle 9 is defined as the 
angle between the momentum vector of the pion from the 
B decay (bachelor pion) and that of the pion of the Dtt 
system in the Dtt rest-frame. 

The x 2 value is calculated using the formula 



"total 



X 



colli 



(7) 



ti 



for cells in a 18 x 18 grid of the two-dimensional his- 
togram. In Eq. (7), ntotai is the total number of cells 
used, N cc in is the number of events in each cell, Nati is 
the expected number of events in that cell as predicted 
by the fit results. The number of degrees of freedom 
(NDF) is calculated as n tota \ — k—l, where k is the num- 
ber of free parameters in the fit. We require A^ t > 10; 
if this requirement is not met then neighboring cells are 
combined until ten events are accumulated. 



C. Matrix element M and Fit Parameters 

This analysis uses an isobar model formulation in 
which the signal decays are described by a coherent sum 
of a number of two-body (Dtt system + bachelor pion) 
amplitudes. The orbital angular momentum between the 
Dtt system and the bachelor pion is denoted here as L. 
The total decay matrix element M for B~ — > D + tt~tt~ 
is given by: 

M = ]T p L e^[N L (x,y)+N L (y,x)} 

L=(0,l,2) 



+ J^Pke*** [A k (x,y)+A k (y,x)], 



(8) 



where the first term represents the S-wave (L = 0), P- 
wave (L = 1) and D-wave (L = 2) nonresonant contribu- 
tions, the second term stands for the resonant contribu- 
tions, the parameters pk and <fik are the magnitudes and 



phases of fc th resonance, while ph and 4>l correspond to 
the magnitudes and phases of the nonresonant contribu- 
tions with angular momentum L. The functions Nl(x, y) 
and Ak(x, y) are the amplitudes for nonresonant and res- 
onant terms, respectively. 

The resonant amplitudes Ak(x,y) are expressed as: 

A k (x,y) = R k (m)F L (p'r')F L (qr)T L (p,q,cos6), (9) 

where Rk(m) is the k th resonance lineshape, F^p'r') and 
Fl (qr) are the Blatt-Weisskopf barrier factors [33] , and 
Tl(jj, q, cos 9) gives the angular distribution. The param- 
eter m (= y/x) is the invariant mass of the Dtt system. 
The parameter p' is the magnitude of the three momen- 
tum of the bachelor pion evaluated in the B-meson rest 
frame. The parameters p and q are the magnitudes of the 
three momenta of the bachelor pion and the pion of the 
Dtt system, both in the Dtt rest frame. The parameters 
p', p, q and 9 arc functions of x and y. 

The nonresonant amplitudes Nl(x, y) with L = 0, 1, 2 
are similar to A k (x, y) but do not contain resonant mass 
terms: 

N (x,y) = 1, (10) 
N^y) = Fi(pV)Fi( g r)Ti( P ,g,cose), (11) 
N 2 (x,y) = F 2 (p'r')F 2 (qr)T 2 (p,q, cos 9). (12) 

The Blatt-Weisskopf barrier factors F^p'r') and 
FL(qr) depend on a single parameter, r' or r, the ra- 
dius of the barrier, which we take to be 1.6 (GeV/c) , 
similarly to Ref. [12]. A discussion of the systematic un- 
certainty associated with the choice of the values of r and 
r' follows below. The forms of Fl(z), where z 
qr, for L = 0, 1, 2 are: 



p'r' or 



Fo(z) = 
Fi(z) = 

F 2 (z) = 



'1 + 4 

1 + z 2 ' 



'9 + 3zg + z 4 
9 + 3z 2 + z A ' 



(13) 
(14) 

(15) 



where zq = p' r' or qor. Here p' and go represent the 
values of p' and q, respectively, when the invariant mass is 
equal to the pole mass of the resonance. For nonresonant 
terms, the fit results are not affected by the choice of 
invariant mass (we use the sum of bid and m v ) used 
for the calculations of p' and go- F° r virtual D* decay, 
D* — > D + tt~ , and virtual B* production in B~ — > B*tt~ , 
we use an exponential form factor in place of the Blatt- 
Weisskopf barrier factor, as discussed in Ref. [12]: 



F(z) = exp(-(z-z% 



(16) 



where z' = rp v for D* — > D + tt~ and z' = r'p v for 
B~ — > B*n~ . Here, we set p v = 0.038 GeV/c, which 
gives the best fit, although any value of p v between 0.015 
and 1.5 GeV/c gives negligible effect on the fitted param- 
eters compared to their statistical errors. 
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The resonance mass term Rk(m) describes the inter- 
mediate resonance. All resonances in this analysis are 
parametrized with relativistic Breit-Wigner functions: 

where the decay width of the resonance depends on m: 



r(m) = r 



^) 2L+1 H Fl( q r), 
q J V m ) 



(18) 



where mo and To are the values of the resonance pole 
mass and decay width, respectively. 

The terms Tl(p, q, cos 8) describe the angular distribu- 
tion of final state particles and are based on the Zemach 
tensor formalism [34]. The definitions of TL(p,q,cos9) 
for L = 0, 1, 2 are: 



T (p, q, cos 9) 
Tx(p 7 q, cos 9) 
T 2 (p,q, cos 9) 



1, 

— 2pq cos 9, 

Ap 2 q 2 {cos 2 9- 1/3). 



The signal function is then given by: 
S{x,y) = \M\ 2 . 



(19) 
(20) 
(21) 

(22) 



In this analysis, the masses of I?* and B* are taken 
from the world averages [5] while their widths are fixed 
at 0.1 McV; the magnitude pk and phase (f>k of the D^ 
amplitude are fixed to 1 and 0, respectively, while the 
masses and widths of Dj resonances and other magni- 
tudes and phases are free parameters to be determined 
in the fit. The effect of varying the masses of D* and B* 
within their errors [5] and widths of D* and B* between 
0.001 and 0.3 MeV is negligible compared to the other 
model-dependent systematic uncertainties given below. 

Since the choice of normalization, phase convention 
and amplitude formalism may not always be the same 
for different experiments, we use fit fractions and relative 
phases instead of amplitudes to allow for a more mean- 
ingful comparison of results. The fit fraction for the fc th 
decay mode is defined as the integral of the resonance de- 
cay amplitudes divided by the coherent matrix element 
squared for the complete Dalitz plot: 

f _ Jpp \Pk{A k {x,y) + A k (y,x))\ 2 dxdy 

h ~ Sm\M\Hxdy ■ (26) 

The fit fraction for nonresonant term with angular mo- 
mentum L has a similar form: 



h 



Jpp \PL(N L (x,y) +N L (y,x))\ 2 dxdy 
J DF \M\ 2 dxdy 



(24) 



The fit fractions do not necessarily add up to unity be- 
cause of interference among the amplitudes. 

To estimate the statistical uncertainties on the fit frac- 
tions, the fit results are randomly modified according to 
the covariance matrix of the fit and the new fractions are 
computed using Eq. (23) or (24). The resulting fit frac- 
tion distribution is fitted with a Gaussian whose width 
gives the error on the given fraction. 



D. Signal Resolution Function 

The detector has finite resolution, thus measured quan- 
tities differ from their true values. For the narrow res- 
onance D 2 with the expected width of about 40 MeV, 
the signal resolution needs to be taken into account. In 
order to obtain the signal resolution on m 2 (Dir) around 
the Z?2 m ass region, we study a sample of MC generated 
B~ — > Xtt~ — > D + n~n~ decays, with the mass and 
width of X set to 2.460 GeV/c 2 (D 2 mass region) and 
MeV, respectively, and subject these events to the same 
analysis reconstruction chain. The reconstructed events 
are then classified into two categories: truth-matched 
(TM) events, where the B and the daughters are correctly 
reconstructed, and self-crossfeed (SCF) events, where one 
or more of the daughters is not correctly associated with 
the generated particle. 

The two-dimensional distribution of cos 9 versus 
m 2 (Dit) for truth-matched events is shown in Fig. 4. 
Since the resolution is independent of cos#, we fit the 
distribution of the quantity q' = m 2 (Dn) — m 2 ruo using a 
sum of two Gaussian functions with a common mean to 
obtain the resolution function for truth-matched events 
(7\Ltm)- The signal resolution for an invariant mass of 
the Dtt combination around the D^ region is about 3 
MeV/c 2 . 

The two-dimensional distribution of cos 9 versus 
m 2 (Dir) for self-crossfeed events is shown in Fig. 5. The 
SCF fraction, /scFj varies from 0.5% to 4.0% with cos6>. 
We fit the /scf distribution with a 4th-order polynomial 
function. The /scf distribution and the result of the 
fit are shown in Fig. 6. The resolution for self-crossfeed 
events varies between 5 MeV/c 2 and 100 MeV/c 2 with 
cos#. We divide the cos 9 interval into 40 bins of equal 
width and use these bins to describe the resolution func- 
tion (T^scf) in terms of a sum of two bifurcated Gaus- 
sian (BGaussian) functions with different means. The 
BGaussian is a Gaussian as a function of q' with three 
parameters, q' the mean, and the two widths, a\ on the 
left and u 2 on the right side of the mean. The form of 
BGaussian is: 



BGaussian(<7' — q' , ai, a 2 ) = 

exp(- 



V2ir(cri+<T2) 
2 

V^Tr(0"l+<T2) 



(q'-Qq) 



if <t < q' ; 



exp( 



(q'-q'o) 



2a. 



¥-) Hq'>q'o, 



(25) 



where q' , 0\ and a 2 are free parameters. 

The signal resolution function is then given by: 

1l{q',cos6) = (1-/scf(cos0)) x1Z TM (q') 

+ /scf(cos0) x K S cF{q',cos9). (26) 

The function lZ{q' , cos 9) represents the probability den- 
sity for an event having the true mass squared m 2 rue to 
be reconstructed at m 2 (Dn) for different cos 9 regions. 

The signal term S in Eq. (4) is convoluted with the 
above resolution function. For each event, the convolu- 
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FIG. 4: Two-dimensional histogram cos# versus m 2 (Dir) of FIG. 5: Two-dimensional histogram cos# versus rn 2 (Dir) of 
the truth-matched events as defined in the text. the self-crossfeed events as denned in the text. 



tion is performed using numerical integration: 
S(x,y)®K= { S{q min + q',q' max ) X 1Z(q' , COS 6) dq' , (27) 



where S is the signal function in Eq. (22), and q min (q max ) 
is the lesser (greater) of x and y. The quantity cos9 is 
determined from q m - ln and g max and is assumed to be 
constant during convolution. The resolution in cos 9 has 
a negligible effect on the fitted parameters. The quan- 
tity q' max is computed using the kinematics of three-body 
decay with q m i n , q' and cos 9. 

The resolution function and the integration method 
in Eq. (27) have been fully tested using 262 MC sam- 
ples with full event reconstruction given below. We have 
compared D invariant mass resolutions for D° — > K~ir + , 
K~ir + ir~ir + and D + — > K~ir + ir + between data and 
MC-simulated events and find that they agree within 
their statistical uncertainties. Estimated biases in the 
fitted parameters due to uncertainties in the signal reso- 
lution function are small and have been included into the 
systematic errors. 



E. Efficiency 

The signal term S defined above is modified in order 
to take into account experimental particle detection and 
event reconstruction efficiency. Since different regions 
of the Dalitz plot correspond to different event topolo- 
gies, the efficiency is not expected to be uniform over the 
Dalitz plot. The term e(x, y) in Eq. (4) is the overall effi- 
ciency for truth-matched and self-crossfeed signal events, 
hence the efficiency for truth-matched signal events is 



(TM(x,y) = e(x,y)(l - /scf(cos0)). 



(28) 



In order to determine the efficiency across the Dalitz 
plot, a sample of simulated B~ — > D + ir~n~ events in 



3? = 



□ 

c/5 
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1 

cos 9 



FIG. 6: /scf(cos#) distribution. The observed self-crossfeed 
fractions (points with statistical errors) are compared to the 
results of the fit (solid curve). 



the Dalitz plot is generated. Some events are gener- 
ated with one or more additional final-state photons to 
account for radiative corrections [35]. As a result, the 
generated Dalitz plot is slightly distorted from the uni- 
form distribution. The number of generated events is 
N g cn =1252k. Each event is subjected to the standard 
reconstruction and selection, described in Section III. In 
addition, we require that the candidate decay is truth 
matched. After correcting for data/MC efficiency dif- 
ferences in particle identification, which are momentum 
dependent and thus vary over the Dalitz plot, the total 
number of accepted events is iV acc = 121,390. We em- 
ploy an unbinncd likelihood method to fit the Dalitz plot 
distributions for generated and accepted event samples. 
The PDF for generated events (PDF gen ) is a fourth-order 
two-dimensional polynomial while the PDF for accepted 
events (PDF acc ) is a seventh-order two-dimensional poly- 
nomial. The efficiency function is then given by: 



, , PDF acc (a;,|/)xJV a 
ctm(^) - x ^ 



(29) 
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FIG. 7: The efficiency for signal decays as a function of 
m 2 (Dn), as determined by MC simulation (points with sta- 
tistical errors) and the results of the fit to the accepted and 
generated distributions (solid curve). 



Fig. 7 shows the efficiency as a function of m 2 (Dir) and 
the fit result for MC-simulated events. 



F. Background 

The background distribution is modeled using MC 
background events, selected with the same criteria ap- 
plied to the data and requiring the B candidate to fall 
into the signal AE region defined in Section III. Events 
in the data AE sidebands could also be used to model the 
background, however in MC studies we find differences 
between the Dalitz plot distributions of the background 
in the signal and sideband regions. Since we find the 
Dalitz plot distributions of sideband events in data and 
in MC simulation to be consistent within their statistics, 
we are confident that the MC simulation can accurately 
represent the background distribution in the signal re- 
gion. Fig. 8a, 8b and Fig. 9a show the Dalitz plot dis- 
tributions of sideband events in data, sideband events in 
the MC sample and background events in the AE sig- 
nal region of the MC sample, respectively. Fig. 8c, 8d 
and 8e show the comparisons of AE sideband events be- 
tween data and MC simulation in ?n mm (-D7r), J72 max (-D7r) 
and m 2 (7T7r) projections, respectively. Here m min (D7r) 
(m 2 nax (£>7r)) is the lesser (greater) of x and y. 

The parameterization used to describe the background 

is: 

B{x,y) = c (^min - <7i) Cl x 

exp (C 2 (<7min - £?l) + C 3 (<?min - <7l) 2 ) 
+ C 4 (g 2 - ?max) C5 X 

exp (c 6 ((j 2 - <7max) + c 7 (g 2 - g max ) 2 ) 
+ c 8 (z - zi)) C£) exp (ci (z - zx) + Cxx(z - zi) 2 ) 
+ ci 5 BGaussian(g max - cx2, ci 3 , cu) 
+ cigBGaussian(z; — ci6, C17, cis), (30) 

where the coefficients Cq to Cig are free parameters to 
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FIG. 8: Comparison of events in the AE sideband: Dalitz 
plot for (a) data and (b) MC-simulated events, and projec- 
tions on (c) m^ nin {D 7r), (d) rn^^^Dit) and (e) m 2 (7T7r) with 
data (points with statistical errors) and MC predictions (his- 
tograms) . 



be determined from the fit, qi = (mj + m^) 2 = 4.04 
GeV 2 /c 4 and q 2 = (m B ~ m^f = 26.41 GeV 2 /c 4 are 
the lower and upper limits of the Dalitz plot, respec- 
tively, zx = (2m,r) 2 = 0.077 GeV 2 /c 4 is the lower limit of 
m 2 (7T7r), <7 m i n is the lesser of x and y, q max is the greater 
of x and y, z is the invariant m 2 (7T7r), and BGaussian is 
given in Eq. (25). 



The projections on m^^Dn), w max (£)7r) and m 2 (7T7r) 
and the result of the fit for the background events in the 
signal region of the MC sample are shown in Fig. 9b, 9c 
and 9d. The x 2 / NDF for the nt is 72/64. 
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FIG. 9: Fit to background events in the AE signal region 
of the MC sample: (a) Dalitz plot and projections on (b) 
m min(Dif), (c) m^^Dn) and (d) m 2 (7T7r) with MC predic- 
tions (points with statistical errors) and the fits (solid curves) . 



RESULTS 



A. Branching Fraction B(B — » D + ir 7v ) 



The total B~ -> D + tt~ 
lated using the relation: 

B = 



7r branching fraction is calcu- 



signal 



[e-B{D+))-2N{B+B-y 



(31) 



where N s { gna \ = 3414 ± 85 is the fitted signal yield given 
in Section III, e is the average efficiency, B(D + ) — (9.22± 
0.21)% is the branching fraction for D + — ► K~tt + it + [5, 
36], and the total number of B + B~ events, N(B + B~) = 
(197.2±3.1) x 10 6 , is determined using N B g and the ratio 
of T(T(4S) -* S+B-)/r(T(4S) -* B°B°) (= 1.065 ± 
0.026) [5]. 

Since the reconstruction efficiencies vary slightly for 
different resonances, the average efficiency is calculated 
by weighing the accepted and generated events by S(x, y) 



with the values for the parameters of our nominal Dalitz 
plot model (discussed below): 



SilT S(xj,yj) x wj 
EST S ( x jiVj) 



(32) 



where Wi is the correction factor which depends on x 
and y due to particle identification efficiency. The value 
e = (8.72 ± 0.05)% is obtained using this method. 

The measured total branching fraction is B(B~ — > 
D+it-tt-) = (1.08 ±0.03) x 10~ 3 , where the stated error 
refers to the statistical uncertainty only. A full discussion 
of the systematic uncertainties follows below. 



B. Dalitz plot analysis results 

The Dalitz plot distribution for data is shown in Fig. 
10. Since the composition of events in the Dalitz plot 
and their distributions are not known a priori, we have 
tried a variety of different assumptions. In particular, we 
test the inclusion of various components, such as the vir- 
tual D* and B* as well as S-, P- and D-wave modeling of 
the nonresonant component, in addition to the expected 
components of D^ , Dq° and background. The D-wave 
nonresonant term does not improve the goodness-of-fit 
and the fraction of D-wave nonresonant contribution is 
close to 0. The results of these tests with variations of 
the models are summarized in Table II. Of these models, 
model 1 produces the best fit quality with the smallest 
number of components, and we choose it as the nom- 
inal fit model. The components considered in this fit 
model are D^ , Dq°, D*, B* and P-wave nonresonant. 
The P-wave nonresonant component is an addition to 
the fit model used in the previous measurement from 
Belle [12]. The sum of the fractions (115 ± 5)% for the 
nominal fit differs from 100% because of destructive in- 
terferences among the amplitudes. The x 2 /NDF for the 
nominal fit is 220/153. To better understand the large 
X 2 /NDF, we look at the contributions to the total \ 2 
from individual cells. We find four cells with \ 2 > 7, 
which inflate the total x 2 - The central points in these 
cells are at: (6.83, -0.722), (6.83, -0.611), (6.83, 0.5) and 



.08, —0.722), where the first value is m^^D-Tr) , and the 
second is cos 8. In order to determine the effect on the 
fitted parameters from these cells, we repeat the nomi- 
nal fit with these cells excluded. The resulting % 2 /NDF 
is 182/149, corresponding to a probability of 3.4%. As- 
suming these large \ 2 contributions are caused by an 
unknown systematic problem, removing them from the 
fit is reasonable. However, under the assumption that 
these high x 2 contributions have a statistical origin, the 
X 2 probability is 0.04% [37]. The low probability indi- 
cates that a model more complex than the isobar model 
may be necessary to describe the characteristics of the 
data. The differences in the fitted and Dq parame- 
ters, when these cells are included or excluded, are as- 
signed to systematic uncertainties, and are much smaller 
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TABLE II: Fit results for the masses, widths, fit fractions and phases from the Dalitz plot analysis of B~ — ► D + -k~it~ 
for different models. The errors are statistical only. The magnitude and phase of the DJ amplitude are fixed to 1 and 0, 
respectively. The background fraction is fixed to 30.4% as described in Section III. The nominal fit corresponds to model 1. 
The labels, S-NR and P-NR, denote the S-wave nonresonant and P-wave nonresonant contributions, respectively. 



Parameter 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


Model 6 


Model 7 


m D ,o (MeV/c^) 


2460.4 ± 1.2 


2460.2 ± 1.0 


2459.1 ± 1.0 


2460.1 ± 1.1 


2461.5 ±1.2 


2458.1 ± 1.1 


2457.4 ± 1.0 


r D »o (MeV) 


41.8 ±2.5 


41.7 ±2.4 


41.1 ± 2.4 


41.8 ±2.4 


42.0 ± 2.5 


41.8 ± 2.4 


41.7 ±2.4 


m n »o (MeV/c 2 ) 


2297 ± 8 


2309 ± 7 


2297 ± 7 


2312 ± 10 


2307 ± 11 


2270 ± 8 


2273 ± 5 




273 ± 12 


285 ±11 


288 ± 12 


289 ± 20 


313 ± 21 


262 ± 1 2 


276 ± 10 




Q9 O 4-1 Q 


0U.0 It l.z 


ol.O It 1.1 


ou. / ± l.D 


39 f\ 4- 1 Q 


*39 QJ.1 Q 

oz.y Zt l.O 


Of) Q 1 1 1 

ou.y it i.i 


4>d»" (rad) 


0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


fo S o (%) 


62.8 ± 2.5 


59.0 ± 2.1 


57.5 ± 1.7 


57.0 ±4.5 


88.0 ± 8.1 


64.8 ± 2.2 


69.7 ± 1.1 


4> D ,o (rad) 


-2.07 ±0.06 


-2.06 ±0.05 


-2.01 ±0.05 


-2.00 ±0.12 


-2.14 ±0.10 


-1.96 ±0.06 


-2.00 ± 0.05 


fn t (%) 


10.1 ± 1.4 


11.3 ±1.5 


9.0 ± 1.2 


11.0 ± 1.5 


9.6 ± 1.3 






4>d* (rad) 


3.00 ±0.12 


2.99 ±0.08 


3.17 ±0.10 


3.05 ±0.12 


2.82 ±0.17 






fB Z (%) 


4.6 ±2.6 


1.4 ±0.5 




1.7 ±0.8 


12.2 ± 5.4 


2.2 ± 1.4 




4>b* (rad) 


2.80 ±0.21 


-2.43 ±0.28 




-2.33 ±0.28 


2.52 ±0.25 


2.28 ±0.38 




/P-NR (%) 


5.4 ±2.4 




1.6 ±0.4 




12.6 ±4.0 


12.7 ±3.1 




^P-NR ( rad ) 


-0.89 ±0.18 




-1.46 ±0.20 




-0.84 ±0.12 


-0.71 ±0.10 




/S-NR (%) 








0.3 ±0.3 


5.2 ±3.8 






^S-NR ( rad ) 








-0.77 ±0.49 


3.30 ± 0.23 






h a (%) 


30.4(fixed) 


30.4(fixed) 


30.4(fixed) 


30.4(fixed) 


30.4(fixed) 


30.4(fixed) 


30.4(fixed) 


NLL 


22970 


22982 


22977 


22982 


22964 


23046 


23125 


x 2 / NDF 


220/153 


240/152 


236/154 


239/153 


216/150 


328/160 


454/161 



than the statistical uncertainties. The removal of these 
cells does not affect the choice of model 1 as the nominal 
fit from Table II. 

Ref. [38] argues for an addition of a Dir S-wave state 
near the Dtt system threshold to the model of the D7T7r 
final state. We have performed tests using the models 
1-4 in Table II with the D* replaced by a Dtt S-wave 
state. Two different parametrizations for Dtt S-wave 
state amplitude are used: one is the function given by 
Eq.(8) of Ref. [38] with the numerator set to constant, 
the other function is the relativistic Breit-Wigner given 
by Eq. (17). Among the tests we have performed with 
these parameterizations, the model with D2, Dq, Dtt S- 
wave (using Eq.(8) of Ref. [38]), B* and P-wave nonres- 
onant gives the best fit with NLL and x 2 /NDF values of 
22997 and 271/151, respectively, which are worse than 
those of the nominal fit even when allowing the Dn S- 
wave's parameters to vary. Each of these models also 
requires large fractions of D^. 

The nominal fit model results in the following branch- 
ing fractions: B(B- -» Dfn~) x B(Df -> D+tt-) = 
(3.5 ± 0.2) x 1CT 4 and B(B- -» D^tt") x B ( d o° -> 
D + ir~) = (6.8±0.3) x 10~ 4 , where the errors are statisti- 
cal only. A full discussion of the systematic uncertainties 
follows below. 

Fig. 11a, lib and 11c show the m^ in (Z)7r), m^ ax (Z)7r) 
and m 2 (7T7r) projections respectively, while Fig. 12a and 
12b show the costf distributions for the Dq° and D^ 3 
mass regions, respectively. The distributions in Fig. 11 
and 12 show good agreement between the data and the 
fit. The angular distribution in the mass region is 



clearly visible and is consistent with the expected D-wave 
distribution of | cos 2 9 — 1/3 | 2 for a spin-2 state. In 
addition, the Dq° signal and the reflection of D^ can 
be easily distinguished in the rn^^Dir) and m^^D-K) 
projection, respectively. The lower edge of m 2 nin (Z?7r) 
is better described with D* v component included than 
without. 

Table III shows the NLL and x 2 / NDF values for the 
nominal fit and for the fits with the broad resonance Dq° 




FIG. 10: Data Dalitz Plot for B~ -> £)+7r"7r". 
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^ 500 1_| i i i i | i i i i | i i i i | i i i i | i i i i | i i i i 




m^(Drc) (GeV 2 /c 4 ) 
«£? 1501 i | i i i | i i i | i i i | i i i | i i i | i i i | 




mLWCGeVV) 



-O 200 1_ i i i | i i i | i i i | i i i | i i i | i i i 




m 2 (ii7i)(GeV 2 /c 4 ) 

FIG. 11: Result of the nominal fit to the data: projections 
on (a) mn lin (D7r), (b) m^ ax (D7r) and (c) m 2 {-K-n:). The points 
with error bars are data, the solid curves represent the nom- 
inal fit. The shaded areas show the DJ contribution, the 
dashed curves show the Dq° signal, the dash-dotted curves 
show the and By signals, and the dotted curves show the 
background. 



excluded or with the J p of the broad resonance replaced 
by other quantum numbers. In all cases, the NLL and 
X 2 /NDF values are significantly worse than that of the 
nominal fit. Fig. 12a illustrates the helicity distributions 
in the Dq° mass region from hypotheses 2-4; clearly the 
nominal fit gives the best description of the data. We 




J -l -0.5 0.5 1 



cos 9 

FIG. 12: Result of the nominal fit to the data: the cosd 
distributions for (a) 4.5 < m 2 (Dir) < 5.5GeV 2 /c 4 region and 
(b) 5.9 < m 2 (DTv) < 6.2 GeV 2 /c 4 region. The points with 
error bars are data, the solid curves represent the nominal fit. 
The dashed, dash-dotted and dotted curves in (a) show the 
fit of hypotheses 2-4 in Table III, respectively. The shaded 
histograms show the cos 9 distributions from AE sidebands 
in data. 



conclude that a broad spin-0 state Dq° is required in the 
fit to the data. The same conclusion is obtained when 
performing the same test on Models 2-5. 



TABLE III: Comparison of the models with different reso- 
nances composition. The labels, S-NR and P-NR, denote the 
S-wave nonresonant and P-wave nonresonant contributions, 
respectively. 



Hypothesis 


Model 


NLL 


X7NDF 




Model 1 (nominal fit) 


22970 


220/153 


1 

2 
3 
4 


D?,D* V ,B*, P-NR 
D* 2 °,D* v ,B:, P-NR, (2+) 
D* 2 °,D* V ,B* V , P-NR, (1-) 
Df,D*,B*, P-NR, S-NR 


23761 
23699 
23427 
23339 


1171/143 
991/144 
638/135 
652/157 
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VI. SYSTEMATIC UNCERTAINTIES 

A. Uncertainties on B(B~ — ► D + iv~ir~) 

As listed in Table IV, the systematic error on the mea- 
surement of the total B~ — ► D + tt~tt~ branching fraction 
is due to the uncertainties on the following quantities: 
the number of B + B~ events in the initial sample, the 
charged track reconstruction and identification efficien- 
cies, and the D + — ► K~-k + ii + branching fraction. The 
uncertainty in the AE background shape, the uncertainty 
in the average efficiency due to the fit models and a pos- 
sible fit bias also contribute to the systematic error. 

TABLE IV: Summary of systematic uncertainties (relative er- 
rors in %) in the measurement of the total branching fraction. 



Systematic Source 


B(B- -^D + tt-tv-) y /0 ) 


Number of B + B~ events 


1.6 


Tracking efficiencies 


2.5 


PID 


1.5 


AE background shape 


1.3 


D + branching fraction 


2.3 


Fit models 


0.7 


Fit bias 


1.0 


Total Systematics 


4.4 



vary the set of cuts on AE, Toes) -R2, cos Oth and mass 
of D + , which increase the number of signal events by 25% 
and the background fraction to 36.5%, and repeat the fits. 
The difference in the fit results is taken as an estimate 
of the systematic uncertainty due to the event selection. 
Fit biases are studied using 1248 parameterized MC sam- 
ples and 262 MC samples with full event reconstruction. 
Small biases are observed for some of the parameters. 
We combine these biases with those coming from high x 2 
cells, as discussed in the previous section, in quadrature 
to obtain the total systematic contribution from the fit 
bias. The uncertainties in PID are obtained by compar- 
ing the nominal fit results with those obtained when the 
PID corrections to the reconstruction efficiency are var- 
ied according to their uncertainties. The uncertainties in 
the efficiency and signal resolution parametrization are 
found to be negligible using the fits to the reconstructed 
MC samples. 

In addition to the above systematic uncertainties, we 
also estimate a mo del- dependent uncertainty that comes 
from the uncertainty in the composition of the signal 
model and the uncertainty in the Blatt-Weisskopf bar- 
rier factors. The model-dependent uncertainties are es- 
timated by comparing the fit results with Models 2-5 in 
Table II and by varying the radius of the barrier, r' and 
r in Eqs. (14)-(16) from to 5 (GeV/c)" 1 . 



The uncertainty on the number of B + B events 
is determined using the uncertainties on T(T(4S) — > 
B+£-)/r(T(4S) -> B°B°) [5] and integrated luminos- 
ity (1.1%). The uncertainty on the input D + branching 
fraction is taken from [36]. The uncertainty in the AE 
background shape is estimated by comparing the signal 
yields between fitting the AE distribution with a lin- 
ear background shape and with higher-order (second and 
third-order) polynomials. The uncertainty in the fit mod- 
els is estimated by comparing the average efficiencies in 
Eq. (32) using Models 2-5 of Table II. The fit bias is es- 
timated to be less than 1% by comparing the generated 
and the fitted value of B(B~ — ► D + tt~tt~) from resonant 
and continuum MC samples. 



B. Uncertainties on Dalitz plot analysis results 

The sources of systematic uncertainties that affect the 
results of the Dalitz Plot analysis are summarized in Ta- 
ble V. These uncertainties are added in quadrature, as 
they are uncorrelated, to obtain the total systematic er- 
ror. 

The uncertainties due to the background parameteri- 
zation are estimated by comparing the results from the 
nominal fit with those obtained when the background 
shape parameters are allowed to float in the fit. The er- 
rors from the uncertainty in the background fraction are 
estimated by comparing the fit results when the back- 
ground fraction is changed by its statistical error. We 



VII. SUMMARY 

In conclusion, we measure the total branching fraction 
of the B~ — ► D + Tr~ir~ decay to be 

B(B~ -» D+tt-j:-) = (1.08 ± 0.03 ± 0.05) x 10~ 3 , 

where the first error is statistical and the second is sys- 
tematic. 

Analysis of the B~ — ► D + tt~h~ Dalitz plot using the 
isobar model confirms the existence of a narrow _D|° and 
a broad Dq° resonance as predicted by Heavy Quark Ef- 
fective Theory. The mass and width of D^ are deter- 
mined to be: 

m D .o = (2460.4 ± f. 2 ± 1.2 ± 1.9) MeV/c 2 and 
r D .o = (41.8 ± 2.5 ± 2.1 ± 2.0) MeV, 

respectively, while for the they are: 

m D ,o = (2297 ± 8 ± 5 ± 19) MeV/c 2 and 
T D ,o = (273 ± 12 ± 17 ± 45) MeV, 

where the first and second errors reflect the statistical 
and systematic uncertainties, respectively, the third one 
is the uncertainty related to the assumed composition 
of signal events and the Blatt-Weisskopf barrier factors. 
The measured masses and widths of both states are con- 
sistent with the world averages [5] and the predictions of 
some theoretical models [39-41]. 
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TABLE V: Summary of systematic uncertainties in the masses, widths and fit fractions of the DJ° and Dq° and the phase of 
Do* . 



Systematic Source 


Am B ,o 


Ar o , 




AT D .o 






A<A n „ 




(MeV/c 2 ) 


(MeV) 


(MeV/c 2 ) 


(MeV) 


(%) 


(%) 


(rad) 


Background parameterization 


1.0 


1.1 


3 


5 


1.2 


0.0 


0.04 


Background fraction 


0.1 


0.4 


2 


1 


0.4 


0.4 


0.00 


Event selection 


0.6 


1.6 


1 


14 


0.3 


0.8 


0.08 


Fit bias 


0.3 


0.7 


4 


8 


0.7 


1.4 


0.02 


PID efficiency 


0.0 


0.1 








0.0 


0.1 


0.01 


Total systematic error 


1.2 


2.1 


5 


17 


1.5 


1.7 


0.09 


Fit models 


1.3 


0.7 


15 


40 


1.5 


17.2 


0.07 


r constant 


1.4 


1.9 


12 


21 


3.8 


7.8 


0.17 


Total model-dependent error 


1.9 


2.0 


19 


45 


4.1 


18.9 


0.18 



We have also obtained exclusive branching fractions 
for Z)?; and Dq° production: 

B(B~ -» Dfn~) x B(Df -» D+tt") 

= (3.5 ±0.2 ±0.2 ±0.4) x 1(T 4 and 

-» DS O x S(D^° -» D+tt") 

= (6.8 ±0.3 ±0.4 ±2.0) x 10~ 4 . 

Our results for the masses, widths and branching frac- 
tions are consistent with but more precise than previous 
measurements performed by Belle [12]. 

The relative phase of the scalar and tensor amplitude 
is measured to be 

(j) D , = -2.07 ± 0.06 ± 0.09 ± 0.18 rad. 

Acknowledgments 

We are grateful for the extraordinary contributions of 
our PEP-II colleagues in achieving the excellent luminos- 



ity and machine conditions that have made this work pos- 
sible. The success of this project also relies critically on 
the expertise and dedication of the computing organiza- 
tions that support BaBAR. The collaborating institutions 
wish to thank SLAC for its support and the kind hospi- 
tality extended to them. This work is supported by the 
US Department of Energy and National Science Foun- 
dation, the Natural Sciences and Engineering Research 
Council (Canada), the Commissariat a l'Energie Atom- 
ique and Institut National de Physique Nucleaire et de 
Physique des Particules (France) , the Bundesministerium 
fur Bildung und Forschung and Deutsche Forschungsgc- 
meinschaft (Germany), the Istituto Nazionale di Fisica 
Nucleare (Italy), the Foundation for Fundamental Re- 
search on Matter (The Netherlands) , the Research Coun- 
cil of Norway, the Ministry of Education and Science of 
the Russian Federation, Ministerio de Education y Cien- 
cia (Spain), and the Science and Technology Facilities 
Council (United Kingdom). Individuals have received 
support from the Marie-Curie IEF program (European 
Union) and the A. P. Sloan Foundation. 



[1] 


N. Isgur and M. B. Wise, Phys. Lett. B 232, 113 (1989). 


[12] 


Belle Collaboration, K. Abe et al, Phys. Rev. D 69, 


[2] 


M. Neubert, Phys. Rept. 245, 259 (1994). 


112002 (2004). 


[3] 


A. F. Falk and M. Luke, Phys. Lett. B 292, 119 (1992). 


[13] 


FOCUS Collaboration, J. M. Link et al, Phys. Lett. B 


[4] 


A. F. Falk and M. E. Peskin, Phys. Rev. D 49, 3320 


586, 11 (2004). 




(1994). 


[14] 


Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 94, 


[5] 


Particle Data Group, C. Amsler et al, Phys. Lett. B 667, 




221805 (2005). 




1 (2008). 


[15] 


CDF Collaboration, A. Abulencia et al, Phys. Rev. D 


[6] 


ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 




73, 051104 (2006). 




232, 398 (1989). 


[16] 


Belle Collaboration, A. Kuzmin et al, Phys. Rev. D 76, 


[7] 


Tagged Photon Spectrometer Collaboration, J. C. Anjos 




012006 (2007). 




et al, Phys. Rev. Lett. 62, 1717 (1989). 


[17] 


ALEPH Collaboration, D. Buskulic et al, Z. Phys. C 73, 


[8] 


CLEO Collaboration, P. Avery et al, Phys. Rev. D 41, 




601 (1997). 




774 (1990). 


[18] 


CLEO Collaboration, A. Anastassov et al, Phys. Rev. 


[9] 


CLEO Collaboration, T. Bergfeld et al, Phys. Lett. B 




Lett. 80, 4127 (1998). 




340, 194 (1994). 


[19] 


DO Collaboration, V. M. Abazov et al, Phys. Rev. Lett. 


[10] 


E687 Collaboration, P. L. Frabetti et al, Phys. Rev. Lett. 


95, 171803 (2005). 




72, 324 (1994). 


[20] 


DELPHI Collaboration, J. Abdallah et al, Eur. Phys. J. 


[11] 


CLEO Collaboration, P. Avery et al, Phys. Lett. B 331, 




C 45, 35, (2006). 




236 (1994). 


[21] 


BABAR Collaboration, B. Aubert et al, Phys. Rev. D 76, 



18 



051101 (2007). 

[22] Belle Collaboration, D. Liventsev et al, Phys. Rev. D 

77, 091503 (2008). 
[23] BABAR Collaboration, B. Aubert et al, arXiv:0808.0333 

[hep-ex], submitted to Phys. Rev. Lett. 
[24] BABAR Collaboration, B. Aubert et al., Phys. Rev. Lett. 

101, 261802 (2008). 
[25] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 

652 (1973). 

[26] Charged conjugate states are implied throughout the pa- 
per. 

[27] GEANT4 Collaboration, S. Agostinelli et al, Nucl. In- 

strum. Methods Phys. Res., Sect. A 506, 250 (2003). 
[28] D. J. Lange, Nucl. Instrum. Methods Phys. Res., Sect. 

A 462, 152 (2001). 
[29] T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 
[30] BABAR Collaboration, B. Aubert et al., Nucl. Instrum. 

Methods A 479, 1 (2002). 
[31] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 

(1978). 

[32] ARGUS Collaboration, H. Albrecht et al, Z. Phys. C 48, 



543 (1990). 

[33] J. Blatt and V. Weisskopf, Theoretical Nuclear Physics, 
p.361, New York: John Wiley & Sons (1952). 

[34] C. Zemach, Phys. Rev. 140, B97 (1965). 

[35] E. Barberio and Z. Was, Comp. Phys. Commun. 79, 291 
(1994). 

[36] CLEO Collaboration, S. Dobbs et al, Phys. Rev. D 76, 
112001 (2007). 

[37] 0.04% is the probability to obtain a \ 2 of 182 or greater 
for a collection of 153 standard normal-distributed num- 
bers, of which the four largest positive deviations have 
been removed. 

[38] D. V. Bugg, arXiv:0901.2217 [hep-ph], to be published in 

J. Phys. G: Nucl. Part. Phys. 
[39] N. Isgur, Phys. Rev. D 57, 4041 (1998). 
[40] D. Ebert, V. O. Galkin and R.N. Faustov, Phys. Rev. D 

57, 5663 (1998). 
[41] Y. S. Kalashnikova and A. V. Nefediev, Phys. Lett. B 

530, 117 (2002). 



